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Abstract. We use Monte-Carlo Markov chain techniques to constrain acceptable parameter regions 
for the Munich L-Galaxies semi-analytic galaxy formation model. Feedback from active galactic 
nuclei (AGN) is required to limit star-formation in the most massive galaxies. However, we show 
that the introduction of tidal stripping of dwarf galaxies as they fall into and merge with their host 
systems can lead to a reduction in the required degree of AGN feedback. In addition, the new model 
correctly reproduces both the metallicity of large galaxies and the fraction of intracluster light. 
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INTRODUCTION 

This paper describes the implementation of a model for dwarf galaxy disruption within 
the semi-analytic framework of De Lucia and Blaizot [1, hereafter DLB07]. The original 
model was suggested by Henriques et al. [2] as a way of both reducing the excess of 
dwarf galaxies and creating the intracluster light (ICL); however this was implemented 
a posteriori, acting only to reduce the dwarf population at the current day. The new 
model follows the stripping of dwarfs as they fall into the halos of their parent galaxy, 
thus gradually reducing their mass: affecting their infall rates, increasing the time-scale 
for, and decreasing the magnitude of, the merger with the central object. 

For the purposes of these conference proceedings, the most important result is that 
the masses of the black holes are reduced, with a corresponding reduction in the level of 
feedback of AGN energy into the interstellar medium. This should be read in conjunction 
with the paper by Chris Short in this volume that investigates the degree to which the 
accretion energy is needed to provide feedback into the intracluster medium (ICM) in 
order to provide the observed entropy excess in clusters. 



METHOD 

Any semi-analytic model has a large number of parameters whose values need to be 
optimized. These are not all independent and correlations between them can give insight 
into the key physical processes that are constrained by the observations. DLB07 has 12 
explicit parameters (as well as lots of hidden ones). We fix all but 6 parameters in our 
analysis: 

• the star formation efficiency, ocsF 



TABLE 1. Statistics from the MCMC parameter estimation for the parameters in the satelHte disruption 
model (Paper 2). The best fit and marginahzed confidence Hmits are compared with the published values 
from DLB07 and with the best fit values obtained without the inclusion of satellite disruption (Paper 1). 
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• the AGN radio mode efficiency, ^agn 

• tlie black liole growtli efficiency, /bh 

• the supernova reheating and ejection efficiency, respectively Cdisk and Chaio 

• the ejected gas reincorporation efficiency, 7ej 

In Henriques et al. [[1, hereafter Paper 1] we introduced a Monte-Carlo Markov Chain 
(MCMC) method for mapping out allowable likelihood regions in parameter space and 
used it to find the best-fitting values for the above parameters, using as constraints the 
^-band luminosity function, the B — V colours, and the black hole-bulge mass ratio. 
Subsequently we have introduced a self-consistent stripping model for the dwarfs [0, 
hereafter Paper 2] . We summarise some of the results from that paper below. 

We build our model on merger trees derived from the Millennium Simulation ([sj], 
using a small but representative subset of the total volume for our MCMC analysis. 
However, the results presented below come from applying the best-fit model to the entire 
simulation. 

RESULTS 

The best fit and confidence limits for the 6 free parameters in the model with satellite 
disruption, together with the published values from DLB07 and the best fit for the model 
without satellite disruption are shown in Table [H The most significant thing to note is 
that, in the disruption model, the star-formation rate is much enhanced and the rate of 
re-incorporation of expelled gas is reduced. This results in a much lower fraction of cold 
gas in the interstellar medium and consequently a lower accretion rate onto the central 
black hole. Thus the black hole masses are smaller in the disruption model, in better 
agreement with observations, and the feedback rates from AGN are reduced. 

Figure [H shows a comparison of the model black-hole/bulge mass relation (contours) 
with observations (crosses). Although the observations are biased to high-mass systems, 
the distribution of systems above and below the best-fit line are correct. 

Although not used to constrain the model, it is gratifying that it produces the correct 
amount of ICL (see Figure 9 of Paper 2). In addition, as shown in Figure |2l it provides 
a much better match to the metallicity distribution of massive galaxies than the original 
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FIGURE 1. The black hole-bulge mass relation for the satellite disruption model (solid contours). The 
red crosses represent observations from Haring and Rix with the best fit to the data points given by the 
red line. 



model (in which massive galaxies accrete too many low-metallicy stars during mergers 
of dwarf satellites.) 

DISCUSSION 

The black hole growth model presented here provides a real challenge for models that 
require AGN heating to raise the entropy of the ICM. Even using the original DLB07 
model, with its higher black hole masses, 35 per cent of the available rest mass energy is 
required (see article by Chris Short in this volume, also Short and Thomas Similar 
conclusions have been reached by Bower et al. |@] using the Durham GALFORM semi- 
analytic model. 
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FIGURE 2. Comparison between the metallicity stars in the best fit for the satellite disruption model 
(solid red lines), in the best fit for DLB07 (dashed red lines) and in observations from Gallazzi et al. |01 
(blue squares and lines). For all the data sets, the central line represents the median value of metallicity in 
each mass bin (the blue squares for the observational data), while the upper and lower lines represent the 
16th and 86th percentiles of the distribution. 
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